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Abstract A technique to find optimal solutions for production-distribution planning in supply chain 
management under intuitionistic fuzzy environment is presented in this paper. The advantages of using 
intuitionistic fuzzy sets over fuzzy sets or crisp sets under imprecise environment are well established. In 1997, 
Angelov proposed optimization technique under intuitionistic fuzzy environment. In 2009, Jimenez and Bilbao 
showed that fuzzy efficient solution may not be Pareto-optimal solution to multiple objective linear 
programming problems in case that one of the fuzzy goals is fully achieved. In 2015, Wu, Liu and Lur redefined 
membership function of fuzzy set theory and proposed another two phase technique. The prime intention to 
maximize the up-gradation of most misfortunate in optimization technique under intuitionistic fuzzy 
environment is better served when some constraints present in existing, well established techniques are 
removed. Moreover, usage of some such constraints present in existing techniques may make a problem 
infeasible. In proposed algorithm new functions are suggested in place of membership functions and non- 
membership functions of intuitionistic fuzzy sets. One standard production-distribution model is taken not only 
to allocate limited available resources and equipment to produce the products over time periods but also to 
determine economic distributors for dispatching product to the retailers. Numerical examples involving this 
production-distribution model further illustrate proposed techniques. Conclusions are drawn at last. 

Keywords Intuitionistic fuzzy optimization • Production-distribution planning • T (+) -characteristic function • T 1 ' 1 - 
characteristic function. 

1 Introduction 

Deterministic optimization techniques are well studied and very popular. But they use crisp set theory and may 
not always represent real life problems. It is very difficult to represent constraints of real life optimization 
problems through crisp constraints. In reality, a small violation of a given constraint is allowed and may 
generate more efficient solution to the problem. But it cannot happen under crisp environment. 

It is known that an optimization problem usually involves objective and constraints. In the last half of century, 
one paradigm shift in mathematics and consequently in optimization techniques is the invention of fuzzy set 
theory by L.A. Zadeh (1965). In optimization problem under fuzzy environment, objective(s), constraints are 
represented by fuzzy sets. Fuzzy optimization techniques are more flexible and generate solutions that are more 
adequate to real life. There exist several techniques in literature to solve optimization problems under fuzzy 
environment. Bellman and Zadeh (1970) introduced decision making in a fuzzy environment. The first 
researchers to publish an operationalization of Bellman and Zadeh’ s proposed approach to 'Decision making in 
a fuzzy environment’ were H. Tanaka and K. Asai (1973) that was published in English in 1974. Zimmermann 
(1976, 1978) further developed mathematical model of optimization problems under fuzzy environment. 
Verdegay (1982) presented another historical paper on fuzzy mathematical programming. Sakawa (1987) 
developed a comprehensive and interactive multiple objective optimization technique under fuzzy environment. 
Recently Jimenez and Bilbao (2009) suggested that fuzzy-efficient solution of multiple objective linear 
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programming problems (MOLPP) may not be Pareto-optimal in case that one of fuzzy goals is fully achieved. 
Their proposed procedure extended two-phase approach of Guu, Wu (1997, 1999) and approach of Dubois, 
Fortemps (1999) to attain Pareto-optimal solutions to MOLPP under fuzzy environment. 

But according to Yan-Kuen Wu et al (2015), proposed approach by Jimenez and Bilbao (2009) cannot 
guarantee to be one general procedure to attain Pareto-optimal solutions. In fact one of the poorly studied 
problems in the field of fuzzy set theory is its membership function itself. Yan-Kuen Wu et al (2015) redefined 
membership functions of fuzzy set theory. But that definition of redefined membership functions and usage in 
mathematical model as well as in given numerical examples in their paper are not analogous. Clearly parts of 
redefined membership functions where values of objective functions do not exceed sum of goals and tolerances 
(for minimization type of objective functions) are used in mathematical model and in given numerical examples. 
But these are not added to set of existing constraints. Thus in problem formulation, redefined membership 
functions are strictly monotonic where as they are not so in definition. Only upper bounds of membership 
functions are removed but there still exist lower bounds of membership functions at zero. And in many cases, if 
the lower bounds of membership functions remain intact, the problem may become infeasible. Membership 
functions of fuzzy sets provide satisficing result when extreme ends of imprecise information can be quantified 
within boundary of zero and one. But it may not be always correct to quantify imprecise information within one 
bounded subset of real line. 

Consider fuzzy set A of all players who may win next US Open grand slam (tennis) tournament. Most of us 
agree with the prediction that Young Donald [current ATP ranking 52] will not be the champion. So it can be 

safely stated that Young Donald is in A with membership value 0. Again chance of Paul Cayre’ s [current ATP 

rank 1000] winning is also nil. So Paul Cayre is also in A with membership value 0. It means that Young 
Donald and Paul Cayre are at same level! But it is certain that even if Young Donald wins, unfortunately Paul 
Cayre probably will not participate and hence will not be the winner. Membership function of fuzzy set theory 
fails to explain this case! An element in fuzzy set may lie partly or never lie or must lie in that set. This concept 
of lying, partly lying or not lying is well measured by membership functions of fuzzy set theory. But it does not 
suit in case of must lying (fuzzy membership values are always one), and never lying (fuzzy membership values 
are always zero). 

On the other hand, fuzzy set theory has been widely developed and several new modifications as well as 
generalizations have appeared. One of these is concept of intuitionistic fuzzy (IF) sets, introduced by K. T. 
Atanassov (1986). Here not only the degrees of acceptance but also degrees of rejection of elements in a set are 
considered so that sum of these values does not exceed unity. Analogous to membership function of fuzzy set 
theory, poorly studied problems in the field of intuitionistic fuzzy set are its membership function and non- 
membership function. Next Plamen Angelov (1997) introduced optimization technique under IF environment by 
using degree of acceptance along with degree of rejection, as a natural extension of optimization technique 
under fuzzy environment. And in literature, the framework of an optimization problem under IF environment is 
usually considered as the following model: 

To maximize the degree(s) of acceptance(s) of intuitionistic fuzzy objective(s) and constraints and to minimize 
the degree(s) of rejection(s) of intuitionistic fuzzy objective(s) and constraints 



Mt7x{/2 ; (x)} ,xe R n ,i = 1,... , p + q (1) 

Min { V- (x)} , x e R n , i = 1 p + q 

x L ' 1 

subject to 

Vj (x)>0,i = l,...,p + q, jU^ix) > v- (x), i = l,...,p + q, 
p- (x) + v- (x) < 1, i = 1, . . . , p + q, x g X. 



Where x denotes decision variables, //■ (x) denotes the degree of membership of x to i th IF set 

Z t (x),i =l,...,p + q , P-(x) denotes the degree of non-membership of x to i 111 IF set z, : (x), i = \,...,p + q , p 
denotes the number of objectives, q denotes the number of constraints, 
X ={xe R" : Ax < b,x>0},b = (b 1 ,b 2 ,...,b m )G R m , A is an m x n matrix and R denotes the set of real 
numbers. Due to symmetry, objective function may be treated as constraint and vice versa. 

The rest of the paper is organized as follows: Sect. 2 introduces definitions of T (+) -characteristic functions and 
T () -characteristic functions . Sect, of 3 highlights loopholes of existing technique and suggests remedy to find 
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optimal solutions to IFO problems. Sect. 4 proposes an interactive and comprehensive algorithm to find optimal 
solutions under IF environment. Sect. 5 introduces a model on production-distribution planning in supply chain 
management under IF environment. Sect. 6 shows through numerical examples how existing technique fails and 
proposed algorithm presents optimal solutions. Sect. 6 also presents two examples where standard membership 
and non-membership functions have not been used. Sect. 7 gives concluding remarks. 

2 Definitions 

It is observed that in case of optimization technique under IF environment, one interesting and useful property 
of membership functions is that higher values of membership functions always give better result of objective 
functions; and analogous property of non-membership functions is that lower values of non-membership 
functions give better result of objective functions. And classical membership functions and non-membership 
functions are not the only candidates having these amazing characteristics. Hence functions viz. T <+> - 
characteristic functions in place of membership functions and ^"’-characteristic functions in place of non- 
membership functions are defined under IF environment. 



Definition 1 

Let S denotes the universal set and A be an IF subset of S. Then T -characteristic function is denoted by 
T^ + \x) and is defined as7L* + * : S — > R that assigns a real number T^} + \x) to each elementxe 5 , here 

T^} + \x) represents degree of membership or acceptance level of x€ S in A . And T^’-characteristic function 

is denoted by T~ / ^ (x) and is defined as T-} ^ : S — > R that assigns a real number T-/ ^ (A) to each element 
A 1 A A' 

( — ) 

X €E S , here T - • (x) represents degree of non-membership or rejection level of x e S in A . 

A 1 



3 Obtaining optimal solution to IFO problems 

IF optimization (IFO) problems are two stages process that include aggregation of objective(s) and constraints 
and then defuzzification. Applying the approach by P. Angelov (1997), analogous to fuzzy optimization 
technique, IFO problem can be transformed into the following problem 

max a- p (2) 

subject to 

p i > a . i = 1 . . . k, v- </?, i = l...k, 
a> p,a + P <\, P>Q, xe X. 

p+q p+q 

where 0C= Min^ix ) , denoting minimal degree of acceptance and p = MaxV i (x ) , denoting maximal degree of 

1=1 i=l 

rejection of IF objective(s) and constraints and k = p + q. It is to be remembered that in problem formulations 
as well as in numerical computations, strictly monotonic parts of definitions of membership functions 
p i (x ), i = 1 , . . . , k and non- membership functions v i (x ), i = \,...,k are used. And triangular IF goal to minimization 

or lesser-equal to type of z, (x) is usually denoted by ( a. ; b. , b '. ) with p. ( z ■ (a.)) = 1 , V. ( z ■ (a.)) = 0, 

p.(z.(b.)) = 0 and V.(z(b .)) = 1 with b. <b. as shown in Fig. 1. Without loss of generality, all constraints are 
assumed to be active here. 

Let us concentrate on philosophy of optimization technique under imprecise environment. In IF optimization 
technique, several operators can be used. Among those, Angelov (1997) has used an operator that is analogous 
to classical max-additive operator. Here prime intention is maximum up-gradation of the most under privileged 
one. Consequently the primary objective in mathematical form is to maximize the difference between minimal 
degree of acceptance and maximal degree of rejection i.e. to maximize a-p . In Angelov’s paper, objective 
function a-p was subjected to the constraints a > p. a + p < p > 0 . But, it is shown below that some of 
these constraints hinder a-p from attaining values higher than unity. Again, as numerical examples shows later 
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Fig. 1 Strictly monotonically extended membership function and non-membership function (not in scale) 



in this paper, some of those constraints may make the problem infeasible. Also, presence of these constraints 
may yield less acceptable optimal solution to the objective function in some cases. These are now discussed one 
after another. 

First consider a simple linear programming problem (LPP) as follows 

max a- P 

subject to 

a > /?, a + /? < 1 , /? > 0 . 

Here a , /? are two decision variables. The optimal solution is a* = 1,/?* = 0. **’ denotes optimality. 

Clearly, maximum value of objective function a- p under given constraints is unity. And a + ft < 1, /? > 0 are 
active constraints in this LPP. Since primary objective is to maximize at-/? under IF environment, therefore in 
order to attain higher value than unity of or-/?, the active constraints of this LPP must have to be removed. 
Again, as numerical example suggests later in this paper, optimal value of z, (jc) cannot be less than <7 ; when 
constraints a + P < 1,/? > Oare present in problem formulation. But, z ; (x) is of minimization or lesser-equal to 
type and IF goals and/or tolerances are imprecise in nature and any value less than to z t {x) is more acceptable 
to decision maker. Hence it is suggested to remove constraints a + P < 1, P > Ofrom model (2). 

Again classical membership functions of IF information have upper bound at unity and classical non- 
membership functions of IF information have lower bound at zero. If constraints CC < 1, p > 0 are added after 
removal of a + P < 1,/? > 0, it does not allow maximum value of at-/? to be higher than unity as well. But, 
primary objective is to maximize a- P under IF environment. Moreover, optimal value of z,(x) cannot be less 
than« ( . But, as discussed earlier, z t (x) is of minimization type and IF goals and/or tolerances are imprecise in 
nature and any value less than a j to z t (x) is more acceptable to decision maker. Hence it is suggested not to add 
constraints OC < 1, /? > 0 to the model. 

On the other hand, constraint or > /? means that minimal degree of acceptance cannot be less than maximum 
degree of rejection. Consequently, Z i (x) cannot attain values greater than C ; ( C ; being abscissa of the point of 

intersection of membership function and non-membership function of Z ; (x)). But as IF goals and/or tolerances 
are imprecise in nature, in many cases, this constraint may make problem infeasible, as shown by numerical 
example. Hence it is suggested to remove the constraint a> /? , if required, so that optimal solution may be 
found. 

Again constraint or > Odoes not allow objective function Z, (x) to take value greater than/?. . But as IF goals 
and/or tolerances are imprecise in nature, it may make the problem infeasible, as shown by numerical example. 
Hence it is suggested to remove constraint a > 0 , if required, so that optimal solution may be found. Hence it is 
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more useful to have simply strictly monotonic functions viz. T (+) -characteristic functions, in place of 
membership functions in IFO technique. 

Again classical definition of non-membership functions of IF set restricts upper bound at unity. If this 
constraint ft < 1 be added, it may make problem infeasible, as shown by numerical example. Therefore it is more 
useful to have strictly monotonic functions viz. ^"’-characteristic functions, in place of non-membership 
functions in IFO technique. 

4 General method to solve IFO problem 

The above ideas can be further integrated into a general framework and an algorithm may be developed to 
obtain optimal solution to optimization problem under IF environment. The steps of the proposed algorithm are 
synthesized as follows 

Step 1. Convert linguistic information to objective data and construct single objective optimization problem 
where objective and / or constraint(s) may be IF in nature. Ask decision maker to specify goal(s) and 
tolerance(s) for IF objective and constraint(s). 

Step 2. Define suitable T (+) -characteristic function in such a way that higher value of T (+) -characteristic 
function always gives better result for IF objective and / or constraint(s). And simultaneously define 
suitable ^"’-characteristic function in such a way that lower value of ^"’-characteristic function gives 
better result for objective and / or constraint(s). It is noted that well-defined T (+) -characteristic functions 
and ^"’-characteristic functions are always finite. 

Step 3. Construct IF optimization problem as 

max a - ft (3) 

subject to 

7j- + (z- (x)) >a,i=l,...XTf ( Zi (x)) <J3, i = 1, . . . , k, 
a>/3,j3> CUe X. 

Step 4. Solve model (3). If it has optimal solution, go to step 8. Otherwise go to step 5. 

Step 5. If model (3) is infeasible, remove constraint a > jB and add constraint a > 0 
becomes 

max a-fS 
subject to 

T ( + (z t (x)) > a, i=l,...,k, T-(zf(x)) < = l,...,k, 

a>0,jS>0,xs X. 

Solve this modified model. If it has optimal solution, go to step 8. Otherwise, go to step 6. 

Step 6. Remove the newly added constraint a > 0 and add the constraint /5 < 1 in last modified model 
modified model becomes 

max a-fi 
subject to 

T i + (z i (x))>a,i=l,...XTr(z i (x))<jB,i=l,...,k, 

f5 < 1, a is unrestricted in sign,.r e X . 

Solve it. If it has optimal solution, go to step 8. Otherwise, remove constraint ft < land solve. If it has 
optimal solution, go to step 8. Otherwise the model is infeasible, even under IF environment; go to step 
7. 

Ask decision maker to modify goal(s) and tolerance(s) for IF objective and constraint(s). Go to step 2. 
It is noted that goalfs) and/or tolerance(s) of IF objectives should better lie within individual maximum 
and minimum values. 

This solution is optimal to single objective optimization problem under IF environment with specified 
goals(s) and tolerance(s). Supply this optimal solution to the decision maker. If he/she is satisfied with 
this solution, stop. 

Otherwise go to Step 7. 



Step 7. 
Step 8. 



. The 
(5) 



. The modified model 

(4) 
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5 Production-distribution planning in supply chain management under IF environment 

With the globalization and evolution of un-nurtured marketplace globally over the time, manufacturing 
companies are forced to design and operate supply chain more efficiently not only to meet customer demands 
but also to maximize the profit under uncertainty. A typical supply chain may consist of a number of suppliers 
and a plant, where processing of the relevant materials adds value to what have been received from suppliers 
and delivers the product to one or more distributors to dispatch them to several distribution centers or retailers. 
And due to variations in soft skills as well as hard skill, transportation cost may vary from one distributor to 
another. Here the aim is to develop a production distribution mathematical model not only to determine the 
production planning but also to select the best contractors. Moreover since maximum manpower available is not 
a crisp number in real life situations, few constraints and objective(s) are IF in nature. The study assumes that 
there is a packaged drinking water bottling factory to produce bottled water. These are then dispatched to 
different distribution centers by several companies as distributors to satisfy demand over a medium planning 
horizon. Also the labor levels may be imprecise in nature. Hence the mathematical model to determine the 
aggregate production plan in a midterm period and to select the economical contracts that dispatch the product 
of supply chain to distribution centers or retailers is as follows 

j : index for destinations j = 1 , ... , J n : index for time periods n = 1 , . . . ,N 

k : index for distributors k = 1 , . . . K Q n :Production volume in period n (boxes) 

I n Tnventory level in period n (boxes) B n :Backorder level in period n (boxes) 

T kii :Number of transported product to destination) by distributor k in period n (boxes) 

1 1 If distributor k dispatches the products to the DCs or retailers 
k [0 Otherwise 

CP n : Production cost per box in period n ($/box) 

CI n Tnventory holding cost per box in period n ($/box) 

CB n : Backordering cos t per box in period n ($/box) 

CT Jkn : Per box transportation cost to destination j in period n by distributor k ($/box) 

FC jhl : Fixed transportation cost to destination) in period n by distributor k ($) 

D n : Demand of product in period n (box) 

R jn : Minimum supply of product for destination) in period n (box) 

MARR: Minimum attractive rate of return as escalating factor for the costs (%) 

U: Required machine-hour to produce one product (machine-hour/box) 

F n : Maximum factory capacity in period n (machine-hour) 

V : Required man-hour to produce one product (man-hour/box) 

M n : Maximum man power available in period n (man-hour) 

S: Required warehouse volume per box of product (m 3 / box) 

AW n : Maximum available warehouse volume in period n (m 3 ) 

Here the objective is to minimize the total cost of the system that includes production, inventory, backorder 
and transportation cost. It is clear that the amounts of product that can be transported in each time period are 
those products that have been produced in the same period plus the amount of products that have been stocked 
in the previous period minus the amount of product that will be maintained in the warehouse for next time 

K J 

period. This leads to the constraint 7 n _, + Q n - I n = XX(n*r j J,Vn.. 

k=l j = 1 

Again it may be assumed that the demand of each distribution center over a special time period might be 
satisfied or backordered, although the backorders have to be fulfilled next period until the end of the planning 
horizon. It leads to the constraint / n l - B n l +Q n - I n + B n = D n , Vn. 

Consequently the amount of product that is produced in each period at least should satisfy the backorder of 
the previous time period. It leads to the constraint Q n > B n l , Vn. 
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And each distribution center for each time period has a minimum requirement of supply that should be 

K 

satisfied. Hence it leads to the constraint X (Y k * T jkn ) > R jn ,\/j,\/n . 

k = 1 

In real life, capacity of man-hour is imprecise in nature. And the usual limitations in the capacity of machine- 
hour, man-hour and the available warehouse volume of the factory are given by the constraints 

U*Q„ < F n Vn, V *<2„ <M n Vn, fj^Y k (S*T jkn )< AW„, Vn. 

*= 1 j=i 

Again the decision maker (DM) can only look for one contractor to dispatch the products to distribution 

K 

centers. Hence the constraint is obtained as X Y k = 1. The non-negativity restrictions on decision variables and 

k = 1 

special restrictions further lead to the constraints Q n ,O n , I n , B n , T jkn > 0 V n and Y k ={0,1}. Here it may be 
assumed that the objective function has imprecise aspiration level. 

Then the mathematical model under IF environment is as follows 



Min 7 = X (CP n * Q n + CI„ * /„ + CB n * S„ + £ X n (CT jh , * T jkn + FC Jhl ))*(! + MARRY 

n= 1 &=1 y'=l 

subject to 

C ln -l n _ l+ Q n -I n =j^(Y k *TjJ,Vn, 



k = 1 7=1 



Q :/ , -B , +Q -I +5 =D , Vn, 

2 n n— 1 n—i n n n 7 7 

: 0. ^ , Vn, C 4 ,„ : X (n * ^ ^ ** . V/', Vn, 

ife=l 

C 5n :U*Q n <F rl , Vn, C 6n :V*Q n <M n , V n, 

C 8 :X F a- = 1 ’ C, : , O n , I n , B n , > 0 Vn, C 10 : Y t ={0,1}. 



C 7 „:XZn(5*^)<AlT„, Vn, 
*=1 7=1 



(6) 



This IF mathematical model is flexible in the value of objective functions and also has the vagueness in some 
constraints. Hence the task is to convert it as model (1) and next as model (2). 

6 Application of production-distribution model under IF environment 

IF mathematical model is now applied on this packaged-drinking water bottling factory. Here the data are taken 
analogous to Ariafar et al (2014). The production planner of the factory has to provide an overall perspective 
regarding the amount of production, inventory and also backordering of their main product for a midterm 
scheduling period, as an aggregate production plan. There are two alternative contractors and three distribution 
centers. The data related to bottling factory are as follows 

• The main product is a bottle of packaged drinking water. Each dozen of bottles is packed in one box. 

• The planning horizon of the factory is three months lune, July and August. 

• The initial inventory is 500 boxes and the end inventory at August is zero. 

• The initial and end back order volume is zero. 

• Production of each box utilizes 0.015 machine hours of the capacity of the factory. 

• The required man power to produce one box of the product equals 0.09 man-hours. The dimension of 

each box is 31x27x37 cm 3 . In 1 m 3 of the space of the warehouse, 32 boxes can be stored. The 
dimension of the warehouse is 15x20 m 2 and height is 1.5m. Hence the volume of the warehouse is 450 
m 3 . 

• The escalating factor of the cost is 8%. 

The other data are given in the following tables. 
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Table 1. Forecasted data on demand and factory capacity 



Table 2. Minimum supply for each 
destination (box) 



Period 


Demand (Box) (D) 


Factory capacity (machine-hour) (F) 


June 


12400 


193 


July 


12000 


185 


August 


12500 


196 



Period 


City A 


City B 


City C 


June 


2000 


2500 


5000 


July 


1500 


2000 


6000 


August 


1000 


2000 


6500 



Table 3. Costs data ($/box) 


Period 


Production cost 


Inventory holding cost 


Back ordering cost 




June 


1.60 


0.10 


0.50 




July 


1.80 


0.08 


0.40 




August 


1.25 


0.09 


0.35 



Table 4. Per box transportation cost for each 
destination ($/box) 



Period 


City A 


City B 


City C 




K=1 


K=2 


K=1 


K=2 


K=1 


K=2 


June 


0.094 


0.112 


0.106 


0.142 


0.064 


0.082 


July 


0.096 


0.118 


0.102 


0.144 


0.060 


0.080 


August 


0.092 


0.112 


0.105 


0.156 


0.058 


0.088 



Table 5. Fixed transportation cost for each destination ($) 



Period City A City B City C 

K=1 K=2 K=1 K=2 K=1 K=2 



June 110 50 180 80 100 30 

July 120 55 175 75 80 25 

August 130 60 160 70 90 35 

Table 6. Imprecise data on maximum manpower available 

in period n (man-hour) 

Period Maximum manpower in period n (man-hour) (M) 



June Availability is 1110 man-hour. It can go up to 1115 man-hour but will not exceed 1118 man-hour. 

July Availability is 1095 man-hour. It can go up to 1100 man-hour but will not exceed 1 104 man-hour. 

August Availability is 1097 man-hour. It can go up to 1100 man-hour but will not exceed 1 102 man-hour. 



Clearly, here J = 3, N = 3 and K= 2. 



6.1. 1. Numerical Example l(removal of constraints (X+[i < 1 ,/J > 0 ) 

In this example, it is shown that the constraints a+j3 < \,(3 > Omay hinder the primary objective CC-fi from 
attaining values higher than unity. 

Suppose that the decision maker chooses goal for the objective function z at $71,000 and decides that cost 
may increase up to $75,000 but in no case it exceeds $76,000. 

Also, the constraints C 61 ,C 62 ,C 63 are IF in nature, as given in Table 6. Based on those information, 
membership and non-membership functions of the IF objective and constraints are constructed as follows 



A(z) = 



1 , 

75000 -z 



4000 

0 , 



if z< 71000 
,if 71000 <z <75000 



0, ifz< 71000 



^i(z) = 



z-71000 

,(/ 

5000 



71000 <z <76000 



1 , 



if z > 75000 



if z > 76000 
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Table 7. Imprecise data on maximum manpower available in period n (man-hour) 

Period Maximum manpower in period n (man-hour) (M) 

June Availability is 1110 man-hour. It can go up to 1115 man-hour but will not exceed 1118 man-hour. 

July Availability is 1095 man-hour. It can go up to 1100 man-hour but will not exceed 1 104 man-hour. 

August Availability is 1077 man-hour. It can go up to 1080 man-hour but will not exceed 1082 man-hour. 



/A (0.09 *0) = 



1, if 0.09*0 <1110 

1115-0.09*0! 

5~ 

0, if 0.09*Q >1115 



,if 111O<O.O9*0<1115 v 2 (0.09 *0,) = 



0, if 0.09*0 <1110 

0.09*0-1110 ^ 1110 < 009 * Gi < 1118 



^(0.09*0,) = 



M 4 (0.09*0) = 



1, 



if 0.09*0, <1095 



1 100 -0 09*0 

^,if 1095 <0.09*0, <1100 v 3 (O.O9*0 2 ) = 

0, 

1, 

1100-0.09 *0 

nr 

0, 



if 0.09*0, >1100 
if 0.09*0, <1097 



if 0.09*0 >1118 
if 0.09*0, <1095 



,if 1095 <0.09*0, <1104 



3 if 1O97<O.O9*03<11OO i/ 4 (O.O9*0 3 ) = 



if O.O9*0 3 >1100 

Then the single objective optimization problem is obtained as follows 



0, 

0.09*0,-1095 . 

9 

X if O.O9*0 2 >1104 

0, if 0.09*0, <1097 

0.09*0,-1097 . 



,if 1O97<O.O9*0 3 <1102 
if 0.09*03 >1102 



max a-j6 
subject to 

ju t >«,Vi = l,2,3,4, v, < /?, Vi = 1,2, 3,4, a+p < 1, a > J3, p > 0, 



(V) 



c,„ : + 0„ - /„ = £ J (X t * T Jkn ), Vn, C 2 „ : - S„_, + 0„ - /„ + B n = D n , Vn, 



k = 1 7=1 



C 3 „:0„>fl_,, Vn, 
Cs„--U*Q„<F n , Vn, 

= l> 






£=1 



C 7 „:ZZn(5*r jfa )<AlT„, Vn, 

/c=l 7=1 

C Q :0 ,0 ,/ ,7\ >0 Vn, 

9 ^-n ’ n ’ n’ 7 «n ’ 



Qo : Xt = {0, 11- 



Using Lingo 15.0.32, the optimal solution is obtained as a*-/3* = 1, z* = $71000. 

As discussed earlier, if constraints a+/3 < > 0 are removed from problem (7), the new optimal solution is 

obtained asef- /? = 1.26,z = $69975.49. Since maximizing up-gradation of most misfortunate i.e. maximizing 
a- p is primary objective, this solution is more preferable to decision maker. Moreover, optimal value of 
objective function is more acceptable after removing constraints a+p < \,p > 0 . Clearly, usage of proposed 
algorithm generates more preferable optimal solution to the problem than traditional IFO technique. 

6.1.2. Numerical Example 2 (removal of upper bound of degree of acceptance) 

After the constraints a+P < \.p > 0 are removed, suggestion to add constraint a <1 may pop up in order to 
restrict upper bound of minimal degree of acceptance at unity. The following example shows that such addition 
is not permissible. 

Suppose that the decision maker chooses goal for the objective function z at $70,800 and decides that cost 
may increase up to $75,000 but in no case it exceeds $78,800. 

Again, constraints C 61 ,C 62 ,C 63 are IF in nature, as given in Table 7. Based on these information, the 
membership and non-membership functions of the IF objective and constraints are constructed. And finally the 
single objective optimization problem is obtained as follows 
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Table 8. Imprecise data on maximum manpower available in period n (man-hour) 

Period Maximum manpower in period n (man-hour) (M) 

June Availability is 1 1 15man-hour. It can go up to 1170 man-hour but will not exceed 1173 man-hour. 

July Availability is 1095 man-hour. It can go up to 1100 man-hour but will not exceed 1110 man-hour. 

August Availability is 1090 man-hour. It can go up to 1100 man-hour but will not exceed 1110 man-hour. 



max a- p 
subject to 

H,>a,Vi = 1,2,3, 4, v,. < /?, Vi = 1,2, 3,4, a < 1, a > p, P > 0, 



( 8 ) 






C, :/ , -B , +Q -I +B =D , Vn, 

2 n n - 1 n - 1 « n n’ ’ 



£=1 7=1 



C 3 „ : Q„ > B n _, , Vn, C 4 : £ (7, * 7^ ) > R jn , V/, Vn, 



/c=l 
K J 



C 5 „ : t/ : * Q„ < F n , Vn, C 7 „ : £ £ Y t (S * T jkn )<AW„ Vn, 

k=l 7=1 

c s =1 ’ C 9 ■Q n ,O n ,I n ,B n ,T jkn > 0 Vn, 



C lo :Y t ={0,l}. 



Using Lingo 15.0.32, the optimal solution is obtained aser*-/?* = l, z* = $70800. As discussed earlier, if 
constraints a < 1, p > 0 are not added, the new optimal solution is obtained as a- P = 1.16, z = $70109.63. 

Since maximizing up-gradation of most misfortunate is prime target i.e. maximizing a- P is primary 
objective, this solution is more preferable to decision maker. Moreover, optimal value of objective function is 
more acceptable when constraints a < 1, P > 0 are not added. Clearly, usage of proposed algorithm generates 
more preferable optimal solution to the problem than traditional IFO technique. 

6.1.3. Numerical Example 3(removal of constraint er > p ) 



This example shows that the constraint a> p may make a problem infeasible whose optimal solution may exist 
with a* < p* . 

Suppose that the decision maker chooses goal for the objective function z at $69,500 and decides that cost 
may increase up to $70,000 but in no case it exceeds $70,300. 

Again, constraints C 61 ,C 62 ,C 63 are IF in nature, as given in Table 8. Based on these information, the 
membership and non-membership functions of the IF objective and constraints are constructed. And finally the 
single objective optimization problem is obtained as follows 

max a- p (9) 



subject to 

ju, > a, Vi = 1,2, 3, 4, v, < p,Vi = 1,2, 3, 4, a > p, 

C ln : + Q n -',,=11 (n * T jkn ), Vn, C 2 „ : - S„_, + Q n - I n + B„ = D n , Vn, 

k=l 7=1 

C 3 „ : Q n > , Vn, C 4 : J (T t * T jkn ) > R jn , V/, Vn, 

k=\ 

C 5n :U*Q„<F n , Vn, C ln HY k (S*T jhl )< AW„ Vn, 

k = 1 7=1 



= 1 , 



C„ :0 ,0 ,/ ,5 ,T. >0 Vn, 

9 ’ n ’ n ’ n ’ tkn ’ 



C 10 :Y,={0,1}. 



Using Lingo 15.0.32, the problem has no feasible solution. 

If constraint a > P is removed, then the optimal solution is obtained as 

a* = 0.14,/?* = 0.53, a* ~P* = -0.39, z* = $69, 928.49. 

Clearly, removal of the constraint a > P on problem formulation and usage of proposed algorithm generates 
optimal solution to the problem where traditional IFO technique fails. 
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Table 9. Imprecise data on maximum manpower available in period n (man-hour) 


Period 


Maximum manpower in period n (man-hour) (M) 




June 

July 

August 


Availability is 1 100 man-hour. It can go up to 1105 man-hour but will not exceed 1 108 man-hour. 
Availability is 1 105 man-hour. It can go up to 1110 man-hour but will not exceed 1111 man-hour. 
Availability is 1077 man-hour. It can go up to 1080 man-hour but will not exceed 1082 man-hour. 


Table 10. Imprecise data on maximum manpower available in period n (man-hour) 


Period 


Maximum manpower in period n (man-hour) (M) 




June 

July 

August 


Availability is 1075 man-hour. It can go up to 1080 man-hour but will not exceed 1083 man-hour. 

Availability is 1 105 man-hour. It can go up to 1110 man-hour but will not exceed 1112 man-hour. 

Availability is 1077 man-hour. It can go up to 1080 man-hour but will not exceed 1082 man-hour. 


6.1.4. 


Numerical Example 4(removal of lower bound of degree of acceptance) 





After the constraint a> J3 has been removed, suggestion to add constraint a > 0 may pop up in order to restrict 
lower bound of minimal degree of acceptance at zero. The following example shows that such addition of 
constraint is not permissible. 

Suppose that the decision maker chooses goal for the objective function z at $69,500 and decides that cost 
may increase up to $70,000 but in no case it exceeds $70,300. 

Again, constraints C 61 ,C 62 ,C 63 are IF in nature, as given in Table 9. Based on these information, the 
membership and non-membership functions of the IF objective and constraints are constructed. And finally the 
single objective optimization problem is obtained as follows 



max a-fd (10) 

subject to 

J u i >ayi = 1,2, 3, 4, v,. <J3 , Vi = 1,2, 3, 4, a> 0, 

K J 

C , : 7 ,+Q -I =YY ( Y*T , ), Vn, C, :/ ,-B ,+Q -1 +B =D , Vn, 

In n-1 n / j / j v k jkn'^ ’ 2n n-1 n-1 *--« n n n’ ’ 

*=1 7=1 

C 3 „ : G„ > B„_, , Vn, C 4y „ : £(1/ * , V/, Vn, 

it=l 

:U*Q n <F n ,V n, C 7 „ : £ ^ n (5 * T jhl )<AW „ , Vn, 

<fc=l 7=1 



C,:£n= 1. 



C a .Q ,0 ,/ ,B ,7. >0 Vn, 

9 ’ n ’ n ’ n ’ 7 KH ’ 



C lo :Y t ={0,l}. 



Using Lingo 15.0.32, the problem has no feasible solution. As discussed earlier, if constraint cr > 0 is not 
added, then the optimal solution is obtained as a* = -0.18,/?* = 0.74, a*-/3* = -0.92, z* = $70,091.42. 

Clearly, proposed algorithm generates optimal solution to the problem whereas restriction on lower bound of 
minimal degree of acceptance makes the problem infeasible. Therefore removal of lower bound of degree of 
acceptance of IF information is proposed. 

6.1.5. Numerical Example 5(removal of upper bound of degree of rejection) 

After the constraint a> f} has been removed and constraint a > 0 is not added, suggestion to add /? < 1 may 
pop up in order to restrict upper bound of maximal degree of rejection at unity. The following example shows 
that such addition of constraint is not permissible. 

Suppose that the decision maker chooses goal for the objective function z at $69,500 and decides that cost 
may increase up to $70,000 but in no case it exceeds $70,300. 

Again, constraints C 61 ,C 62 ,C 63 are IF in nature, as given in Table 10. Based on these information, the 
membership and non-membership functions of the IF objective and constraints are constructed. And finally the 
single objective optimization problem is obtained as follows 










12 



Table 11. Imprecise data on maximum manpower available in period n (man-hour) 



Period 


Maximum manpower in period n (man-hour) (M) 


June 

July 

August 


Availability is 1075 man-hour. It can go up to 1080 man-hour but will not exceed 1083 man-hour. 
Availability is 1 100 man-hour. It can go up to 1105 man-hour but will not exceed 1 107 man-hour. 
Availability is 1072 man-hour. It can go up to 1080 man-hour but will not exceed 1082 man-hour. 



max a-fd 
subject to 

//. > a,\fi = 1,2, 3, 4, v. < j3,Vi = 1,2, 3, 4, jd< 1, 



(ID 



C, :I ,+Q -1 =yy (Y. * T . ), Vn, C, :7 ,-B ,+0 -7 +B =D , Vn, 

In n-1 ^~-n n / j / j v k jkn'^ 7 2n n-1 n-1 n n n’ 7 



£=1 y=l 



C in -.Q„>B n _ t , Vn, 
C 5n :U*Q„<F n , Vn, 

C 8 :£n=l> 



C,jn-t( Y k * T jJ S7? ; „,Vy,Vn, 



fc=l 
AT / 



Ci n :y?,Y k (S* T j J<AW n , Vn, 

*=1 7=1 

C q :Q ,0 ,/ ,£ ,r.. >0 Vn, 

y ^n 7 n 7 n 7 n 7 7ati 7 



Cio : Xt = (0,1}- 



Using Lingo 15.0.32, the problem has no feasible solution. As discussed earlier, if constraint /5 < 1 is not 
added, then the optimal solution is obtained as a * = -0.75, /3* = 1.09 ,a*-/3* = -1.84,z* = $70,371.21. 

Clearly, proposed algorithm generates optimal solution to the problem whereas restriction on upper bound of 
maximal degree of rejection makes the problem infeasible. Therefore removal of upper bound of degree of 
rejection of IF information is proposed. 



6.2.1. Numerical example using proposed algorithm with first set of I <+> -characteristic functions and T* 1 - 
cliaracteristic functions 



Step 1. Suppose that the decision maker chooses goal for the objective function z at $69,500 and decides that 
cost may increase up to $70,000 but in no case it exceeds $70,300. Again, constraints C 61 ,C 62 ,C 63 are IF in 
nature, as given in Table 11. 

Step 2. Based on these information, T (+) -characteristic functions and T^-characteristic functions are constructed 
as follows 



r(+) 70000-z 

500 

1080 -0.09*0, 

r (+) mna*/o i — ^V0 3 



r +, 2 (O.O9*0 1 ) = 



r <+> 3 (0.09*0,) = 



1105 -0.09* 0, 



V0 2 



z-69500 

800 

, , 0.09*0-1075 

r'-', (0.09*0) = s V0 



r <_) 3 (0.09*0,) = 



8 

0.09*0,-1100 



V0 2 



7’ <+) 4 (O.O9*0 3 ): 



1080-0.09*03 



V0 3 



7% (0.09* ft): 



0.09*0,-1072 

10 



V0 3 



Step 3. Then the single objective optimization problem under IF environment is obtained as follows 
max a-jd 



subject to 

T i + (z i (x))>a,i= 1,2,3, 4, 7jV( z .(jc)) < J3,i = 1,2, 3, 4, a > p,p > 0, X. 

Step 4. Using Lingo 15.0.32, the problem has no feasible solution. 

Step 5. Next remove constraint a > /3 and add constraint ctr > 0 and solve the modified model. Using Lingo 
15.0.32, this problem also has no feasible solution. 
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Step 6. Next remove constraint a > 0 and add constraint P < 1 . Solve this modified model by using Lingo 

15.0. 32. It has no feasible solution. As discussed earlier, remove constraint J3 < 1 and solve. Using Lingo 

15.0. 32, the optimal solutionis obtained as «* = -(!. 69,/?* = l. 21, a*-/?* = -1.90,z* = $70,345. 83. . 

Step 8. This solution is optimal to the problem under IF environment for given goals and tolerances. Supply this 
solution to the decision maker. If he/she is satisfied with this solution, stop. Otherwise, request the decision 
maker to choose different goals and tolerances for IF objective and constraints. 

Hence proposed algorithm generates optimal solution to the problem under imprecise environment where as 
traditional IFO technique fails. 

6.2.2. Numerical Example using proposed algorithm with second set of I i+> -characteristic functions and l i ~ > - 
characteristic functions 

Suppose decision maker is unable to choose goals and tolerances for imprecise objective and constraints. In that 
case, strictly monotonically decreasing T (+) -characteristic functions and strictly monotonically increasing T 1 " 1 - 
characteristic functions can be used for minimization and lesser equal to type of objectives and constraints. 
Consider the following example. 

Step 1. Suppose that the decision maker is unable to choose goals and tolerances for imprecise objective z and 
constraints C 61 , C 62 , C 63 . 

Step 2. Construct strictly monotonically decreasing T <+) -characteristic functions and strictly monotonically 
increasing ^"’-characteristic functions of IF objective and constraints as follows 

T i+ \(z) =-z Vz , r l_) (z) = z Vz, T< +, 2 (0.09*Q) = ■ -0.09*Q VQ, 

r (_) 2 (0.09 *Q) = 0.09 VQ, (0.09 *e 2 )= -0.09 *e 2 Vg 2 , r <_, 3 (0.09 *e 2 ) = 0.09*2, V2 2 , 

r (+> 4 (0.09*2 3 ) =-0.09*2, V2 3 , r H 4 (0.09*2 3 )=0.09*2 3 V2 3 . 

Step 3. Then the single objective optimization problem under IF environment is obtained as follows 
max a- p 

subject to 

7j + (z-(x)) >a,\= 1, 2, 3. 4, Tf (z : f (x)) < p , i = 1 ,2, 3, 4, a > p, P > 0, x e X . 

Step 4. Using Lingo 15.0.32, the optimal solution is obtained as z* = $68,565.85. 

Step 8. This is the optimal solution to the IFO problem. 

Here completely new set of T (+) -characteristic functions and ^"’-characteristic functions have been used for 
imprecise objective and constraints. Several other such T <+ ’-characteristic functions and ^"’-characteristic 
functions may be constructed. These can produce optimal solutions to the IFO problems. 

7 Conclusion 

In this paper, an algorithm is developed to find optimal solutions for optimization problems under IF 
environment. IF technique is one of the richest apparatus for formulation of optimization problems under 
imprecise environment, thereby generating more satisficing result than crisp optimization technique. 
Philosophically the prime objective of optimization problem is maximum up-gradation of most under-privileged 
imprecise information, i.e. in mathematical form the objective is to maximize the difference between minimal 
degree of acceptance and maximal degree of rejection. 

In this paper it is shown how constraints present in traditional well established technique hinders objective 
function from attaining values higher than unity. Again, constraint that minimal degree of acceptance must not 
be less than maximal degree of rejection may make a problem infeasible etc. 

Therefore it is proposed to use strictly monotonic T (+) -characteristic function in place of membership function 
and strictly monotonic T^’-characteristic function in place of non-membership function in IFO technique. This 
algorithm is designed in such a way that it always generates optimal solution under IF environment, whenever 
crisp feasible space is non-null. 

In this paper, proposed algorithm is applied on a production distribution planning in supply chain 
management. The information supplied being imprecise in nature, IF optimization technique is used to find 
optimal solution. The proposed algorithm generates optimal solution to this model whereas traditional IFO 
technique fails for given goals and tolerances. Moreover it is interactive in nature. Hence scope of further 
improvement is always present in it. 
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It may further be concluded that, under these circumstances, the issue of getting no feasible solution under IF 
environment yielded from several published methods seems worthwhile or even necessary to reconsider. 

Moreover, a decision maker may not be present always to supply goal(s) and/or tolerance(s). He/she may not 
be well equipped as well. Moreover, it may be sometimes difficult to construct classical membership and/or 
non-membership functions. Usage of strictly monotonic T <+) -characteristic functions and strictly monotonic T < ‘ > - 
characteristic functions helps in these cases also. Many such T (+) -characteristic functions and ^‘’-characteristic 
functions can be constructed as necessary to generate optimal solution under imprecise environment. 
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